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The optical conductivity of graphene and bilayer graphene in quantizing magnetic fields is studied. 
Both dynamical conductivities, longitudinal and Hall's, are analytically evaluated. The conductiv- 
ity peaks are explained in terms of electron transitions. Correspondences between the transition 
frequencies and the magneto-optical features are established using the theoretical results. The main 
optical transitions obey the selection rule with An = 1 for the Landau number n. The Faraday 
rotation and light transmission in the quantizing magnetic fields are calculated. The effects of 
temperatures and magnetic fields on the chemical potential are considered. 
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I. INTRODUCTION 

The most accurate investigation of the band structure 
of metals and semiconductors is a study of the Landau 
levels through experiments such as magneto-optics^ - — 
and magneto-transporiji^— . In magnetic fields, the clas- 
sical and quantum Hall effects are observed, as well as 
the polarization rotation for transmitted (Faraday's ro- 
tation) or reflected (Kerr's rotation) lights. However, the 
interpretation of the experimental results involves a sig- 
nificant degree of uncertainty, because it is not clear how 
the resonances can be identified and which electron tran- 
sitions they correspond to. 

Comprehensive literature on the graphene family can 
be described in terms of the Dirac gapless fermions. Ac- 
cording to this picture, in graphene, there are two bands 
at the K hexagon vertexes of the Brillouin zone with- 
out any gap between them, and the electron dispersion 
can be considered as linear in the wide wave-vector re- 
gion. For the dispersion linearity, this region should be 
small compared with the size of the Brillouin zone, i.e. 
less than 10 8 cm , providing the small carrier concen- 
tration no <C 10 16 cm~ 2 . Pristine graphene at zero tem- 
perature has no carriers, and the Fermi level should di- 
vide the conduction and valence bands. However, un- 
doped graphene cannot be really obtained, and so far 
purest graphene contains about uq ~ 10 9 cm -2 of car- 
riers. Then the following problem appears — how do 
Coulomb electron-electron interactions renormalize the 
linear dispersion and does graphene become an insulator 
with a gap? 

Semiconductors with the gap are needed for electronic 
applications. Investigations of the graphene bilayer and 
multilayer are very popular as the gap appears when the 
bias is applied. Here Slonczewski, Weiss, and McClure 
(SWMC) should be mentioned because they have stated 
the description of a layered matter 1 — with interactions 
strong in a layer and weak between layers. The theory 
contains several parameters which are the hopping in- 
tegrals for nearest neighbors. Such a picture has been 
examined in many experiments^. 



The theoretical solution for the band problem in mag- 
netic fields often cannot be exactly found. A typical 
example is presented by graphene layers. For bilayer 
graphene and graphite, the effective Hamiltonian is a 
4x4 matrix giving four energy bands. The trigonal warp- 
ing described by the small parameter 73 in the effective 
Hamiltonian provides an evident effect. Another impor- 
tant parameter is the gate-tunable bandgap U in bilayer 
graphene. In this situation, the quantization problem 
cannot be solved within a rigorous method. To over- 
come this difficulty, several methods have been proposed 
for approximate 8 -^ - —, numerical^ 2 - - — , and semiclassical 
quantization^ - —. 

In this paper, our attention is focused on the dynamic 
conductivity of monolayer and bilayer graphene in the 
presence of a constant magnetic field in z— direction. We 
consider the collisionless limit when the electron collision 
rate is much less then the frequency of the electric field. 
Then, the accurate theoretical results can be obtained for 
Faraday's rotation and transmittance through graphene 
layers. The present paper is organized as follows. In Sec. 
II we recall the electron dispersion in the monolayer and 
bilayer graphene. In Sec. Ill we describe in detail the 
quantization in magnetic fields. In Sec. IV the longi- 
tudinal and Hall conductivities as well as the Faraday 
rotation are described. Effects of temperatures and mag- 
netic fields on the chemical potential are considered in 
Sec. V. Section VI contains a summary of the discussed 
results. 



II. ELECTRON DISPERSION IN MONOLAYER 
AND BILAYER GRAPHENE 

Electron dispersion in graphene. 

The symmetry of the K point is C$ v with the threefold 
axis and reflection planes. This group has twofold rep- 
resentation with the basis functions transforming each 
in other under reflections and obtaining the factors 
exp(±27rz/3) in rotations. The linear momentum dis- 
placements from the K point, taken as p± = ^fipx —Py, 
transform in a similar way. The effective Hamiltonian is 
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invariant under the group transformations, and we have 
the unique possibility to construct the invariant Hamil- 
tonian linear in the momentum as 



H(P) 



vp + 
vp_ 



(1) 



where v is a constant of the velocity units. The 
same Hamiltonian can be written using the tight-binding 
model. 

The eigenvalues of this matrix give two bands 



£1.2 = TV 



Pi 



Tvp, 



where the subscript s = 1,2 numerates these two bands 
(holes and electrons). The gapless linear spectrum arises 
as a consequence of the symmetry, and the Fermi energy 
coincides with the band crossing (the Dirac point) due to 
the carbon valence. The cyclotron mass has the form 



m(e) 



1 dS(e) 
2tt de 



and the carrier concentration at zero temperature 
n(ep) = £ 2 F /ith 2 v 2 is simply expressed in terms of the 
Fermi energy ep. 

Tuning the gate voltage, the linearity of the spectrum 
has been examined in the Schubnikov-de Haas studies 3 ^ 
with the help of the connection between the effective 
mass and the carrier concentration at the Fermi level 
m(ef)u = ^Ti^J~Kn(eF)- The "constant" parameter v 
was found to be no longer constant. At low carrier con- 
centrations n ~ 10 9 cm~ 2 , it exceeds its constant value 
V = 1.05 ± 0.1 x 10 8 cm/s for concentrations n > 10 11 
enr 2 by the factor of 3. 

This is a result of electron-electron interactions which 
become stronger at low carrier concentrations. The log- 
arithmic renormalization of the velocity was found by 
Abrikosov and Beneslavsky^I for the three-dimensional 
case and in Refs . 32 ' 33 for two-dimensional graphene. No- 
tice, that no phase transition was revealed even at low- 
est carrier concentration. We have also to conclude that 
the Coulomb interactions do not create any gap in the 
graphen spectrum. 

We recall the peculiarity of graphene conductivity in 
the absence of the magnetic fiel d 34 ' 35 . For the optical 
frequency range, when the spacial dispersion of conduc- 
tivity is not significant, the intraband electron transitions 
make a contribution 



2ie 2 T . . , a . 

— In (2 cosh il), (2) 



■kTi(lo + it 1 ) 



which has the Drudc-Boltzmann form at the large chem- 
ical potential /i T. 

At the zero temperature, the interband electron con- 
tribution can be presented in the simple form 



v ' Ah 



( s (u + 2m) 2 

(w - 2/i) - — In — - 

2n [w — 2li) z 



where the 9— function expresses the threshold behavior 
of interband electron transitions at u) = 2/i. The temper- 
ature smooths out all the singularities in this formula. In 
high frequency region u 3> (T,li), the interband transi- 
tions make the leading contribution into conductivity 



ef 
4ft' 



having the universal character independent of any mate- 
rial parameters. This frequency region is limited above 
by the band width of around 3 eV. 

Making use the universal conductivity, one can calcu- 
late the light transmission through graphen o 36 ' 37 in the 
approximation linear in conductivity 

47r e 2 
T = 1 Re (j(lu) cos 9 = 1 — ir— cos 9 , 

c he 

where 9 is the incidence angle of light. In excellent agree- 
ment with the theory, for the wide optical range, several 
experimental groups^ 8 - - — observe the light transmission 
through graphene as well as bilayer graphene where the 
difference from unity is twice as larger. It is exception- 
ally intriguing that the light transmission involves the 
fine structure constant a = e 2 /he of quantum electrody- 
namics having really no relations to the graphene physics. 

For the frequenciy range, where the intraband term 
plays the main role, the plasmon excitations are 
possibl o 34 ' 41 with the dispersion 



— ^ln(2cosh^). 



and relatively small damping, determined by the electron 
relaxation r _1 . The plasmon has the same dispersion, 
k 1 / 2 , as the normal 2d plasmon. However, it shows the 
temperature dependence at low carrier concentrations, 
fj, < 2T. 

Electron dispersion in bilayer graphene. 

Bilayer graphene has attracted much interest partly 
due to the opening of a tunable gap in its electronic 
spectrum with an external electrostatic field. Such a phe- 
nomenon was predicted in Refs i 42 ' 43 and was observed in 
optical studies controlled by applying a gate voltage 4 ^—. 

The Hamiltonian of the SWMC theory can be 
writte n 24 ' 25 near the K points in the Brillouin zone in 
the form 



H(p) = 



U 
vp- 
7i 



vp+ 
U 



Hvp+flo 73WP-/70 



71 74VP-/70 
74VP-/70 73WP+/70 
— U vp^ 
vp + —U 



(3) 

where p± = ^fip x — p y . The nearest-neighbor hopping 
integral 70 ~ 3 eV corresponds with the velocity param- 
eter v = 1.5ao7o = 10 6 m/s and the in-layer inter-atomic 
distance a = 1.415 A . The parameters 73,4 ~ 0.1 eV de- 
scribe the interlayer interaction at the distance <io=3.35A 
between layers. 
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Hamiltonian (|3]) give four levels labeled by the number 
s = 1,2,3,4 from the bottom. For {7 = 0, the twofold 
degeneration £2 = £3 exists at p x — p y = 0, as a conse- 
quence of axial symmetry. The parameter U is included 
in the bilayer Hamiltonian to describe the gate voltage. 
At U y^z 0, the gap appears between the £2 and £3, and 
these bands acquire the form of " mexican hat" . 

Two vertexes, K and K', in the Brillouin zone are 
transforming each in other under reflection. Such the 
reflection changes the U sign. Therefore, the levels at 
these vertexes do not coincide in the presence of a gate. 
The levels at K' point can be obtained from levels at the 
K point by changing their signs. 



III. GRAPHENE IN MAGNETIC FIELDS 

In the presence of the magnetic field B, the momen- 
tum projections p+ and p- become the operators with 
the commutation rule {p+,p-} = —2ehB/c. We use the 
relations 



vp+ — ujb a, vp- = a 



+ 



involving the creation a + and annihilati on a ope ra- 
tors with the energy parameter ljb = vy/2\e\hB/c = 
36.2^/B [ Tesla] meV . 

For graphene, one seeks the eigenfunction of Hamilto- 
nian ([1]) in the form 



CO*) 



(4) 



where ip n (x) are orthonormal Hermitian polynomials 
with the numbers n > 0, including the exponential factor 
with one of the space coordinates in the Landau gauge. 
The number s = 1,2 numerates the solutions at given 
n. We write only one of two x, y space coordinates. The 
degeneracy proportional to the magnetic field will be in- 
cluded in the final results. Every row of the matrix (p} 
turns out proportional to the definite Hermitian polyno- 
mial which can be canceled from the eigenvalue equa- 
tions. We obtain a system of linear equations 



-e WflV" 



r 1 

^ sn 

a 2 







(5) 



for the eigenvector C sn with two eigenvalues, s = 1, 2, 



e sn = ^u B \/n 

at given n = 1, 2... 

The wave function columns write 



ct _ 1 / 1 



V2{ 



and 



(6) 



(7) 



for s = 1, 2, correspondingly. If n — 0, there is only one 
level £10 = with C\ a — 0, Cf = 1 as follows from Eqs. 
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FIG. 1: Landau levels in doped bilayer graphene with the 
Fermi energy ef — 120 meV in the magnetic field B = 10 T 
and the gate voltage U = ±100 meV at the K,K' points of 
the Brillouin points (right and left panels, respectively); the 
arrows show the electron transitions from the 1 23) level, where 
n = 3 is the Landau number and s =2 is the level number at 
given n from the bottom. 



For bilayer graphene, we seek the eigenfunction of 
Hamiltonian (J3|) as a column 



C 1 



(8) 



similar to Eq. . The every row in Hamiltonian ([3]) be- 
comes again proportional to the definite Hermitian func- 
tion, if the terms with the trigonal warping 73 are omit- 
ted. These terms can be considered within the perturba- 
tion theory or the semiclassical approximation. 

Canceling the Hermitian functions from the equations, 
we obtain a system of the linear equations for the eigen- 
vector C s „ 

/ 




V 0^4 \fri — 1 



where the band number s — 1,2,3,4 numerates the so- 
lutions at given n from the bottom, u>b — vW2\e\hB/c 
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FIG. 2: Longitudinal (a xx , the left panel) and Hall (ff xv , 
the right panel) dynamical conductivities in the real and 
imaginary parts for doped graphene with the Fermi energy 
e F = 120 meV in the magnetic field B = 10 T. 



and C04 = 74Wb/7o- 

The eigenvalues of the matrix in Eq. (J9)), see Fig. [T] 
are easily found using the personal computer. Without 
74, we have for e sn the simplified equation 



[(U-e sn ) 2 ~u 2 B n}[(U+s sn ) 2 -Lo 2 B (n-V}+j!(U 2 



0. 



For each Landau number n > 2, there are four eigen- 
values e sn and four corresponding eigenvectors, Eq. ([5]l. 
marked by the band subscript s; we will also use the 
notation \sn) for levels. In addition, there are four lev- 
els. One of them exists for n = with the eigenvector 
C = (0,1,0,0) as is evident from Eqs. flS), ©. Using 
the perturbation theory, we take the trigonal warping 
into account 17 and obtain 



ei(n = 0) = U + 



( WB73 V 

V 70 / 



|C?< 3 I 



U-e s >(3) 



(10) 



Other three levels indicated with n = 1 and s = 1, 2, 3 
are determined by the first three equations of the system 
© with Cf r = 0. The |21) level is very close to the |10) 
level. 

Alternatively, the semiclassical quantization can be ap- 
plied for relatively weak magnetic fields when the cy- 
clotron frequency ^ff- is small compared to the Fermi 



energy. Then, we get the Bohr-Zommerfeld condition in 
the form 



ehB 



S(e) = 2tt 



T 

4 



6(e) 



(11) 



Here S(e) is the cross-section area of the classical electron 
orbit in the p x ,Py space for the energy e; n is an integer 
supposed to be large. The integer T is the number of 
the smooth turning points on the electron orbit. For 
instance, there are two smooth turning points for the 
Landau quadratic model and only one for the skipping 
electrons reflected by the hard edge. 

If spin is neglected, 5 = and T = 2 for the Landau 
quadratic model, and 5 = 1/2 and T = 2 for monolayer 
graphene. In these two cases, the semiclassical result 
coincides with the rigorous quantization. 

Notice, that the r5(e)-phase depends on the energy and 
it is closely connected with the topological Berry phase 52 . 
The (5-phase was evaluated for bismuth in Ref^, and it 
was considered again for bismuth in Ref^. For graphite, 
semiclassical quantization was applied in Refill. How- 
ever, in the general case, the evaluation of the 5— phase 
still attracts a widespread interest - — . 

In the simplest case of bilayer graphene without trigo- 
nal warping, we find the Berry phased 



5(e) 



-eU 



-eU 



q 2 — e 2 



U 2 



y/4U 2 e 2 + (e 2 - C/ 2 ) 7l 2 



For the ungaped bilayer, [7 = 0, the Berry phase 6(e) = 
0. The Berry phase depends on the energy and 5 = =Fl/2 
at e = ±[7, correspondingly. At the larger energy, e ^> U, 
the Berry phase 5 —> ^fU/ji. The effect of the trigonal 
warping is considered in Ref.— . 



IV. MAGNETO-OPTICS IN GRAPHENE 
LAYERS 

An important peculiarity of conductivity in magnetic 
fields is an appearance of the Hall component a xy (oj). 
The Hall conductivity violates the rotation symmetry of 
graphene around the major z— axis. This implies the ro- 
tation of the linear polarized electromagnetic wave, i. e., 
the Faraday and Kerr effects for transmitted and reflected 
waves, correspondingly. 

First of all, the electron transitions are possible be- 
tween the levels with the neighboring Landau numbers n 
and various bands s, and therefore the resonance denom- 
inators A ss ' n = e sn — e s ', n +i arise in the conductivity 
tensor. 

Calculations 17 give the conductivities for layer 
graphene in the collisionless limit, when the electron col- 
lision frequency V is much less than the level splitting, 



,(u) 

'xy(u) 
X [(lu + lT - 



ia x 



44 



w £i X — 'J S3 II I j 

l(J 0— T" lC-7-\ d 



TV 



2 ^ Ass'n 

n,s,s' 

1 ± (lu + iT - A s 



(12) 
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FIG. 3: Contribution of K point in dynamical conductivities 
for doped bilayer graphene with the Fermi energy ef = 120 
meV in the magnetic field B = 10 T; the interlayer hoping 
integral 71 =360 meV. 



Here A/ ss /„ = /(e a / n +i) - f{e sn ) is the difference of the 
Fermi functions f{e sn ) = [exp( £ "^T M ) + and 

ri , — c 2 c 1 -i- r 13 r 14 

"ss'n — ^sn^s'n+l ^ ^sn^s'n+l 
+ (C S 1 nC s 'n+l + CsnCs' n+l)74/70 

is the dipole matrix element expressed in terms of wave 
functions C s „, Eqs. ([7]), ©, and a — e 2 /4fi is the 
graphene universal conductivity. 

The electron transitions obey the the selection rule 

An = 1 . 

Besides, the renormalization of the dipole moments 
due to trigonal warping results in weak lines with the 
selection rule 

An = 2. 

Then, we get an additional term in Eq. (|12p by substi- 
tuting 

d ss 'n = (73/7o)Cs„Cs'„+2 

as the matrix element and replacing the subscript n+1 — > 
n + 2. We have to notice, that the 74 corrections give 



FIG. 4: Contribution of K' point in dynamical conductivities 
for doped bilayer graphene; the parameters are the same as 
in Fig. El 

the linear (in the small parameter 74/70) contribution to 
the conductivities at the main electron transitions with 
An = 1. The 73 corrections are quadratic, however, they 
result in an appearance of new weak resonant transitions 
with An = 2. 

For graphene, the calculated conductivities are shown 
in Fig. [2]in the frequency range 0+0.5 eV at the magnetic 
field B =10 T. The electron parameters are taken as 
follows: the Fermi energy Ep =120 meV and the electron 
scattering rate V =5 meV. The peaks in absorption, the 
left panel in Fig. [21 correspond to the electron transitions 
between the levels. The strongest line at 113 meV is 
exited by the 1 10} — > |21) transitions. Other lines are 
doublet exited by transitions of the type |ln) — > \2,n+l) 
and |l,n + 1) — > |2,rt) for n from 1 to 4. All lines obey 
the selection rule An = 1. 

For bilayer graphene, the system of lines is more com- 
plicated. At K point (the gate voltage U =100 meV), it is 
shown in Fig. Q] schematically and in Table for the Fermi 
energy ep — 120 meV. The contributions of K and K' 
points into conductivities are compared in Figs. [3] and [4] 

The absorption resonances in Fig. [31 the left panel, are 
composed by the following electron transitions: 

16 meV: |35) -+ |36), 

254 meV (very weak):|25) -+ |36), 

285 meV: |26) -+ |37), |27) -+ |36), 
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TABLE I: Landau levels \sn) for 2 < n <12 in bilayer; other 



levels 


are 


10) = 


100, |11>= 


= -384, |21) 


= 82, |31) 


= 401 in meV. 


\ n 
s \ 
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10 
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-124 
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-155 
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-186 
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113 
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FIG. 5: Transmittance (T, the left panel) and Faraday rota- 
tion angle (Qf, the right panel) for doped graphene with the 
Fermi energy ef — 120 meV in the magnetic field B — 10 T. 



301 meV: 1 10) - 
317 meV: |27) -. 
348 meV: |28) - 

142), 
380 meV: |210) - 
398 meV (weak) 
410 meV: \2ll) 

|46),|34)-> 
439 meV 
468 meV 
496 meV 



131), 

|38),|28 >^ |37), 

|3Q>, |29> |38>, |33> — ► |42), |21) 



|39),|29)^|310),|32)^|43), 
|33) -► |44), 

-> |310), |210> -> |311),|35) 
|45),|35>^ |46), 
|212) ->■ |311), |211) -» |312), 
|213) ->■ |312), |212) |313), 
|214) ->■ |313), |213) |314), |22) |31). 
Graphene and bilayer graphene affect £/ie transmission 
and the Faraday rotation to the linear order in the fine 
structure constant a as well as the reflected light in- 
tensity quadratic in a. The conductivities a xx (to) and 
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FIG. 6: Transmittance (T, the left panel) and Faraday rota- 
tion angle (Qf, the right panel) for doped bilayer with the 
Fermi energy ef = 120 meV in the magnetic field B — 10 T. 



o~xy(w) allow calculating the Faraday rotation and the 
transmittance as functions of the frequency. Because the 
conductivity of the layers are small, we can use the linear 
approximation in a. The transmission coefficient T and 
the Faraday angle Of for the free standing layers write 
as 



47T 27T 

T = 1 Rca xx (u) , 9f = — Rea xy (uj). (13) 

c c 



Results of calculations are shown in Fig. [5]for graphene 
and in Fig. [S] for bilayer. We take into account that 
the points K and K' have the different electron levels 
in bilayer at U ^ 0, and both these points contribute 
independently in observable quantities. 

It is evident that the interpretation of the Faraday ro- 
tation governed by the conductivity a xy (ui) is much more 
complicated in comparison with the transmittance con- 
trolled by the longitudinal conductivity. 

The positions of the lines for fields in the range of 10 
- 30 T agree with observations of Refsii^. 
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V. EFFECT OF TEMPERATURES AND 
MAGNETIC FIELDS ON THE CHEMICAL 
POTENTIAL 

In the previous figures, we assume that the tempera- 
ture and level width V are much less than the level split- 
ting |e S)n +i — £ s ,n |- As is known, the chemical potential 
fi in semiconductors changes while the temperature in- 
creases. We know also that de Haas- van Alphen oscilla- 
tions can be observed at low temperatures. We can find 
the amplitude of these effects assuming that the carrier 
concentration Uq should be constant when the tempera- 
ture or the magnetic field rises. 

For simplicity reasons, let us consider now graphene, 
where e sn = Tw B V", n = 0, 1, 2.., s = 1,2, g = 4, 
dn/de s = 2e s /ui 2 B , u>b — v \j2eh B/c. Electrons in the 
upper band and holes in the lower band influence each 
other, especially at high temperatures. In order to in- 
corporate electrons (fi > 0) and holes (ju < 0) in a sin- 
gle scheme, let us introduce the variable e„ = ±e sn for 
s = 2, 1 correspondingly. Thus the carrier concentrations 
for 2d systems write 



n = sign(^) 



geB 



2ithc ^— ' 

71=0 



[f(e n - n) - /(£„ + /i)] 



where f(s n — fi) is the Fermi-Dirac function, g is the 
spin and valley factor, and no is considered as positive 
whereas fi can be of both signs. 

When the level splitting is far less than the chemical 
potential, we can use the Poison summation formula 
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FIG. 7: Chemical potential versus temperature for the carrier 
concentration 10 10 cm -2 ; the exact solution to Eq. (jTSJ) in 
solid line, asymptotes for low, Eq. (|17[) . and high, Eq. (|18l) . 
temperatures (dashed and dashed-dotted lines, correspond- 
ingly). 



We find the temperature dependence of the chemical 
potential differentiating the equation (|15l) with respect 
to the temperature: 







ede 



dy{e) dfj, dy{e) 



dfi dT dT 



(16) 



n 



2irhc 



poo pOO 

/ dny(e n ) + £ / dne 2 ™ kn y(e r 

Jn=0 k ^ Jn=0 



(14) 

with y(e n ) = \f{e n — fi) — f(e n + The second term 
in brackets goes to zero at B — > and then the first term 
gives the carrier concentration without the magnetic field 



At low temperatures, T \fi\, we can evaluate these 
integrals. Because the integrand in the first integral is 
proportional to the Dirac S— function, it gives fi^jtp- The 
second integral equals 7r 2 T/3. Then, we get the equation 

dfj, 



7T 2 



o 



no(£^0) = 



geB 
2nhc 



dedn / 'dey(e) 



This term depends on the magnetic field only through 
the chemical potential: 



n (B -> 0) = 



!J 



2irv 2 h 2 Jo 



edey(e) . 



(15) 



At given carrier concentration no determined by dop- 
ing or a bias voltage, the equation (Ti"5j) gives the temper- 
ature dependence of the chemical potential (see Fig. [7]) 
and, particularly, the Fermi energy at zero temperature, 
M (T = 0) = e F : 



e F = ^%v\fA-Kn Q /g, 
where =p stand for holes and electrons correspondingly. 



with a solution 



= ei - —T 



or at low temperatures 



fi = e F - —T 2 /e F . 
6 



(17) 



At high temperatures, T> the integral, Eq. (|15l) . 
gives the approximate dependence 



I Ml 



Trno(hv) 2 
41n2T 



(18) 



shown in Fig. [7] While the temperature grows, the 
chemical potential tends to its value /i = in undoped 
graphene. 
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The second term in brackets (|T4|) is easily calculated 
at low temperatures: 



/>oo 

1(B) = £ / deee 2 ^ kn y(e) 



^ ' Ql 



sin(27rfcn M ) 



(19) 



k>0 



sinh(27r 2 fc<T) : 



where n M = n 2 /uj 2 b is the Landau number, corresponding 
to the Fermi energy and n^ = 2u/oj^ is its derivative. 
Owing to the denominator, this term is relatively small 
at weak magnetic fields, 47r 2 |u|T ^> w|, and we accept 
this condition. 

Then, we consider the effect of the magnetic field on 
the chemical potential in the way similar to Eq. (| 16|) : 



= 



dB 



ede 



dy(e) , dl(B) 
9/i 



dB 



(20) 



The integral here equals Integrating the obtained 
differential equation with the initial condition n = e F at 
B = T = and taking Eq. (fl~7|) into account, we get the 
dependence of the chemical potential on the magnetic 
field and the temperature 



VI. SUMMARY AND CONCLUSIONS 



In conclusions, we discuss the dynamic conductivity 
of monolayer and bilayer graphene in the optical range, 
where the ac frequency is much larger than the electron 
scattering rate. The trigonal warping in bilayer can be 
considered within the perturbation theory at strong mag- 
netic fields larger than 1 T approximately. For weak mag- 
netic fields, when the Fermi energy much larger than the 
cyclotron frequency, the semiclassical quantization with 
the Berry phase included can be applied. The main elec- 
tron transitions obey the selection rule An = 1 for the 
Landau number n, however An = 2 transitions due to the 
trigonal warping with the small probability are also pos- 
sible. The SWMC parameters are used in the fit taking 
their values from the previous dHvA measurements. The 
calculated conductivities, longitudinal and Hall's, permit 
to evaluate transmittance and Faraday's rotation of light 
in the graphene layers. The agreement between the cal- 
culations and the measured Faraday rotation and trans- 
missions in graphene in the quantizing magnetic fields 
is achieved. Assuming that the carrier concentrations is 
fixed by gate voltage or dopant, we find the effect of tem- 
peratures and magnetic fields on the chemical potential. 



/i = e F - I(B)/e F - — T 2 /e F , 
6 



(21) 



where we have to take e F instead of /i in 1(B) (|19p for 
the case of the relatively large \e F \ ujb- 

The equation (|2"Tj) represents de Haas-van Alphen os- 
cillations and the chemical potential shift in 2d systems. 
We emphasize that the considered effects have no rela- 
tions to the electron-electron or electron-phonon inter- 
actions but come out purely from the broadening of the 
distribution function. The corresponding impact, pro- 
portional to the temperature, on the width T is consid- 
ered by Ozerin in Ref^. 
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